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Numerical calculation of vector electromagnetic modes of plano-concave mi- 
crocavities reveals that the polarization-dependent reflectivity of a flat Bragg 
mirror can lead to unexpected cavity field distributions for nominally paraxial 
modes. Even in a rotationally symmetric resonator, certain pairs of orbital an- 
gular momenta are necessarily mixed in an excitation-independent way to form 
doublets. A characteristic mixing angle is identified, which even in the paraxial 
limit can be designed to have large values. This correction to Gaussian theory 
is zeroth-order in deviations from paraxiality. We discuss the resulting nonuni- 
form polarization fields. Observation will require small cavities with sufficiently 
high Q. Possible applications are proposed. 

Unconventional beam profiles in free-space paraxial optics have recently received renewed 
attention, based in large part on the new degrees of freedom made accessible in light-matter 
interactions when Gaussian beams carry not only spin but orbital angular momentum 
[1, 2]. To realize the full potential of the resulting rich phenomenology in quantum optical 
applications [3] , it is an important next step to turn from free-space optics to ultrahigh-Q 
resonators. In this Letter, we study numerically a type of cavity that enables robust and 
spectrally addressable creation of nearly paraxial light fields with orbital and polarization 
profiles that are, surprisingly, not predicted by the standard solutions of paraxial theory. 
We build upon the initial examination in [4, 5]. 
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Beams carrying orbital angular momentum require an axially symmetric environment, 
and hence the cavity in this study has a rotation axis z. The dome-shaped resonator 
is shown in Fig. 1 (a). Standard paraxial wave solutions with definite orbital angular 
momentum are the Laguerre- Gauss beams (LG). They are labeled by the number p of 
radial nodes, and the orbital angular momentum quantum number I: in polar coordinates 



where L p (x) is an associated Laguerre polynomial and w is the beam waist radius. All 
beams of the same order N = 2p + \£\ form a degenerate manifold in the sense that their 
longitudinal wave functions are the same for a given frequency u. The degree of degeneracy 
is g — 2(iV + l), including two possible polarization states (e.g. right and left circular, with 
Jones vectors a ± ); this allows the formation of linear combinations of Laguerre-Gaussians 
to obtain a wide variety of transverse beam profiles. The cross-sectional electric field within 
a degenerate LG manifold of order N can be expanded as 



where I = —N + 2j (j = 0, 1, . . .N), and p is fixed by 2p + \£\ = N. This carries 
over to paraxial resonators where the discrete resonant frequencies u Vj n are labeled by a 
longitudinal index v and the mode order N, and do not depend on £ and p individually. 

Our exact numerical solutions reveal that corrections [6] to the standard paraxial res- 
onator theory stated above lead to a splitting of the iV-multiplets in (2). Compared to the 
equidistant transverse mode spacing in u v> n (governed by the Guoy phase shift), this addi- 
tional non-paraxial splitting typically occurs on a scale so fine that it can be disregarded in 
large resonators. However, ultrahigh-finesse microcavities, which are now becoming tech- 
nologically feasible, will make it necessary to quantify this substructure. The problem is 
then analogous to zeroth-order degenerate perturbation theory, the small quantity being 
the paraxiality parameter, h = X/(ttw) = ta.n(divergence angle). 

The question how the degeneracy in (2) is lifted goes beyond paraxial theory. A first 
guess would be that the new eigenmodes are the LG basis modes LG e p a s where s = ±1 is 
the spin. All modes are then uniquely labeled by p, I and s. This is indeed what we find 
when modeling the planar mirror in Fig. 1 ( perfect electric or magnetic conductor 

(-^tangential = Or -^tangential = 0). 

In this Letter, however, we focus attention on the case where the planar mirror is a realis- 
tic Bragg stack. Then, most of the dome cavity modes are not labeled by a unique (p, £, s) , 
even for the smallest values of the paraxiality parameter h. What makes distributed Bragg 
reflectors (DBR) different is their polarization-dependent plane- wave reflectivity r s / p (6), 
where s/p denotes linear polarization perpendicular/parallel to the plane of incidence and 
9 is the angle of incidence. Figure 1 (b) illustrates that reflection phase shifts 4> s / p with 
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Figure 1: (a) Schematic cross section of the cavity. The dome is metal, the base is either 
metal or Bragg mirror. The radius of curvature R is ten times larger than the 
length L. (b) Polarization-dependent reflectivity phases S / P = arg(r a / p ) of the 
different mirrors used in the calculation, versus angle of incidence 9. (c) Index 
profiles of the quarter- wave Bragg mirrors used in Figs. 2 and 3. Each has 36 
pairs of layers of index 3.0 and 3.5. 



strong 6 dependence may occur even if the DBR design wavelength A<j is chosen so that 
the modes of interest are at the center of the stop band, where |?"s/p(#)| ~ 1- We have 
calculated the fully vectorial cavity modes, using a recently developed numerical technique 
combining a basis function expansion in the dome region with the transfer matrix tech- 
nique for the DBR [4, 5]. What we describe below can be called DBR-induced, paraxial 
spin-orbit coupling of light. 

The numerical method is not restricted to near-paraxial modesfl], and the results pre- 
sented here are observed over a wide range of h; but we shall focus on the limit h ^ 1 
where {LG e p a s } should constitute a suitable basis in which to express the (transverse) 
cavity fields. The vectorial modes can be chosen to be eigenf unctions of the total angular 
momentum around the z axis, with integer eigenvalues m. If orbital and spin angular 
momenta i and s are "good quantum numbers" , then m = i + s. As illustrated in Table 1 
of Ref. [4], specifying N and m in (2) singles out pairs of LGpO" 5 , with t = m ± 1 (unless 
\m\ = N + 1, in which case only \£\ = N occurs). We shall call these pairs doublets because 
their paraxial degeneracy is in fact lifted. This is to be distinguished from a remaining 
reflection-induced exact degeneracy between modes differing only in the sign of m. By 
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Figure 2: Mixing angle a for numerically calculated vector cavity modes, versus (a) relative 
detuning between the mean wavelength A of the doublet C, D and the Bragg 
design wavelength Ad, and (b) S — P . Line fits use the function given in the 
text. Cavity dimensions are L = lO/xm and R = 100/mi. 



fixing m, the latter degeneracy can be disregarded for our purposes. 

The lowest order for which the above doublets exist is N = 2, and we consider this 
case from now on. Both A = LG'j'o" 4 " and B = LGgcr - are transverse basis functions with 
m = l. If m is a good quantum number but I and s are not, then the transverse parts of the 
actual m = 1, N — 2 cavity modes (denoted symbolically by C, D) will be superpositions 
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(3) 



where a is a mixing angle: when a = 0, the modes are well approximated by pure LG 
profiles. Among the different mirrors shown in Fig. 1, only the conductor shows a = in 
general. The doublet's resonant frequencies satisfy u d — uj c — > for h — > 0, but the mixing 
angle a generally does not vanish in this paraxial limit. 

Figure 2 (a) gives the variation of a as Ad is changed. The sigmoid shape of the curves 
means that in a typical cavity, one will find the doublet C, D at a value of \a\ near 45°. 
This is the furthest from pure LG modes we can get, in spite of the rotational symmetry of 
the whole cavity. Note that a = 90° again describes pure LG states, and in going through 
a total variation by Act = 90° the modes exchange character. This dramatic change occurs 
over a narrow A — Ad interval, going through the "unmixed" state a — at which point 
the resonance frequencies uq, wd exhibit an avoided crossing (not shown). 

The non-conservation of orbital and spin angular momentum does not rely on paraxiality: 
the x-y field components of a cavity eigenstate with total angular momentum m are linear 



4 



Figure 3: Cross-sectional electric field of C and D for the point marked by the arrow in 
Fig. 2 (a) (mirror design 2, A = 401. lnm, divergence angle 9 m 3.6°). The 
ellipses represent only the x-y components of the electric field; white or black 
filling indicates sense of rotation. Each polarization ellipse is decorated with a 
line pointing along the instantaneous field vector to indicate the phase. Here the 
mixing angle is a = 39°, and uj-q — ujq = 2ir x 3.4GHz. Inset: E x in the x-z plane. 

combinations of two opposing circular polarization fields which individually transform un- 
der coordinate rotations according to different orbital angular momenta i — m ± 1; modes 
A and B in (3) are the paraxial realization of this general fact. A Bessel-wave expansion[ I] 
of the mode shows that one of the orbital wave functions can in fact be made to vanish 
exactly if r s (6) = r p (6) for the planar mirror. This occurs for pure magnetic or electric 
conductors. In the DBR cavity, on the other hand, we only have r s = r p at 9 = 0, cf. Fig. 
1 (b). Since even the most paraxial modes contain nonzero 8 in their angular spectrum, 
one cannot generally factor the field into a fixed polarization vector multiplying a scalar 
orbital wave function. Within the conventional paraxial approach, where this factorization 
is performed at the outset, the consequences of a finite a at nearly vanishing h are lost. 

By changing the detuning of the doublet from Ad, we also change the phase difference 
S — P at any finite angle of incidence 9. As discussed elsewhere [7], one can identify an 
"effective" 6 ~ h at which to evaluate S — P . Plotting a versus this variable S — P in 
Fig. 2 (b), the universality of the spin-orbit coupling is revealed: the data differ widely 
in the mirror designs and wavelengths used, but collapse onto very similar curves. The 
broad appearance of the dotted curve in 2 (a) arises mainly because the reflection phases 
of mirror 3 are less sensitive to detuning. The data are fit by the functional], tana = 
— (0s — P — 0o)/(r + |0s — P — 0o|)- The offset 0o for attaining circular polarization 



(a = 0) accounts for the fact that 4> s and 4> p are not equal at zero detuning when evaluated 
at finite 9. The crossover between a ~ ±45° has width T, and persists in the zeroth-order 
limit h — > where r — > 0. 

In Figure 3, we illustrate the polarization patterns produced by the spin-orbit coupling. 
Polarization ellipticity and orientation vary spatially. In particular, the single mode C 
exhibits circular polarization near the axis but radial linear polarization within its second 
radial maximum; in mode D, a crossover from circular to azimuthal linear polarization 
occurs. The L = 10/im cavity height in this example allows uj-q — uoq to be resolved despite 
the finite widths of the modes, which is taken into account in our numerics (cf. the DBR 
design in Fig. 1, and the data on decay widths in Ref. [4]). Assuming that both mirrors 
have a power reflectivity r 2 in the range (1 — r 2 ) < i o _3 - 5±0 - 5 ; we estimate that the effects 
shown here should be observable for cavities with paraxiality 0.09 < h < 0.2 if the height 
L of the resonator lies in the window 12A < L < 60A. The mode patterns in Fig. 3 are 
governed by the mixing angle a; Fig. 2 suggests the intriguing possibility of externally 
controlling a by changing the cavity dimensions or Bragg wavelength Ad- 

When addressed by laser pulses at u>c and u^, the doublet modes could interact with 
one, two, or more asymmetric (polarization-sensitive) quantum dots embedded in the DBR. 
Thus the mixed modes may allow a new scheme of quantum processing. Cavities of the 
type studied here can also act as filters or laser emitters of pure or mixed LG beams, useful 
as "optical tweezers" [2] or as carriers of information [8]. 

Provided that at least one of the mirrors is a Bragg stack, any stable cavity should 
exhibit paraxial mixing of orbital angular momenta. Its observability hinges on the ability 
to distinguish true degeneracies from "quasi" -degeneracies caused by the breakdown of 
Eq. (2). Our quantitative analysis shows that the necessary requirements can be met by 
realistic cavities that are being pursued for quantum-optics applications. 

We thank Mike Raymer for valuable discussions. This work was supported by NSF 
Grant ECS-0239332. 
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